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We propose and investigate a novel, minimal, and experimentally testable framework for baryo- 
genesis, dubbed dexiogenesis, using baryon number violating effective interactions of right-handed 
Majorana neutrinos responsible for the seesaw mechanism. The distinct LHC signature of our 
framework is same-sign top quark final states, possibly originating from displaced vertices. The 
region of parameters relevant for LHC phenomenology can also yield concomitant signals in nucleon 
decay experiments. We provide a simple ultraviolet origin for our effective operators, by adding a 
color-triplet scalar, which could ultimately arise from a grand unified theory. 


Gauge singlet right-handed neutrinos (RHNs) provide 
perhaps the simplest explanation of non-zero neutrino 
masses, as demanded by a large and well-established 
body of neutrino oscillation data, thereby allowing the 
Standard Model (SM) to be a renormalizable theory of 
Nature. If neutrinos are Dirac particles, the associated 
RHNs only couple to the SM via Yukawa interactions of 
negligible strength, yjg ;$ 10 —12 , and are not expected 
to be directly detectable in the foreseeable future. On 
the other hand, if the observed neutrinos are Majorana 
states, they most naturally get their mass from a see¬ 
saw mechanism pQ. In that case, RHNs may have (D( 1) 
couplings to SM neutrinos, but they would then have to 
be exceedingly heavy, > 10 14 GeV, far beyond the reach 
of terrestrial experiments. However, it may very well be 
that RHNs are much lighter, and they can be directly 
probed in lepton number violating processes in collider 
experiments E3U and perhaps other searches, such as 
those for proton decay [7j. 

In this letter, we assume that the RHNs associated 
with the seesaw mechanism are near the weak scale 
< 1 TeV. Higher dimensional operators involving the 
RHNs are generically present. In particular, we will fur¬ 
ther assume that the RHNs have baryon number vio¬ 
lating interactions, mediated by dimension 6 operators 
involving right-handed quarks and suppressed by a scale 
~ 1 — 10 TeV. We will show that these assumptions al¬ 
low for direct generation of a baryon number asymmetry 
through RHN decays in the early Universe, which we dub 
dexiogenesis (dexios: Greek for the right hand). This is 
in contrast to canonical leptogenesis [SJ where the lep¬ 
ton asymmetry needs to be further processed into baryon 
number through electroweak sphalerons [9]. Our direct 
baryogenesis mechanism is most constrained by nucleon 
decay bounds. However, we show that for viable parame¬ 
ters one could have distinct collider signatures. This sce¬ 
nario, with dim-6 operators, is an effective field theory 
and can be embedded in a simple renormalizable model, 
and possibly a grand unified theory (GUT), where addi¬ 
tional signals are expected to arise at the Large Hadron 


Collider (LHC) and future colliders. For a partial list of 
other works whose subjects have some overlap with that 
of this letter, see, for example, Refs. [HHI7J. 

Our point of departure is the SM augmented by two 
Majorana RHNs N a , a = 1, 2, the minimum required for 
a realistic seesaw mechanism based on current neutrino 
data. We add the following terms to the Lagrangian 

C N = M a N c a N a + y%HN a Li + h.c., (1) 

where M a is the mass of N a , i = 1,2,3 enumerates SM 
generations, and y ai is a Yukawa matrix; H and U are 
the Higgs and lepton doublets of the SM, respectively. 

Light neutrino masses m„ <0.1 eV, implied by the 
oscillation data, can be generated from the renormaliz¬ 
able interactions in Eq. Q, via the seesaw mechanism: 
(m ,,)ij ~ yf}y^ (H) 2 /M a . Nonetheless, the SM, hence¬ 
forth defined to include the Lagrangian in Eq. 0 , is 
widely expected to be an effective theory that is further 
enriched with new interactions at higher scales. This ex¬ 
pectation is strongly motivated by the need for a dark 
matter candidate and also a baryogenesis mechanism to 
generate the observed baryon asymmetry of the Universe. 

Baryogenesis requires a source of baryon number viola¬ 
tion JTS] . Hence, it may be necessary to extend the SM by 
effective operators that violate baryon number jl]3]. Such 
operators are suppressed by scales associated with new 
physics, often considered to be very high, > 10 15 GeV, 
as implied by nucleon decay constraints. However, it is 
compelling to look for scenarios where new physics arises 
at lower scales. For one thing, it is reasonable to assume 
that the physics underlying the Higgs potential is not 
very far from the weak scale. Such physics would then 
have the added benefit of being potentially testable. 

Motivated by the above considerations, we assume the 
following baryon-number violating operators involving 
the RHNs, in addition to those made up of only the ob¬ 
served SM fields, 

y^ijk i^ilm 

^BV = ~K 2 -[N aU idjdk)R + -^ 2 -[N a di\ii[QiQm\L + h.c. , 

( 2 ) 
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grams in Fig. [T] will lead to a baryon asymmetry s = 
r(7Vi —> tbb) — r(iVi —> tbb )/( 2 rjy-j), where the width of 
N\ is given by 


= 


lAipMf 

1024tt 3 A 4 


T K/Mf) 


( 4 ) 


FIG. 1: Tree and two-loop diagrams for dexiogenesis. 


where i,j,k are family numbers of right-handed quark 
(u, d) mass eigenstates and l,m enumerate left-handed 
quark ( Q ) generations. Here, X l J k and K^ m are generally 
complex constants determined by the ultraviolet (UV) 
theory. These operators could arise from grand unified 
theories, as shown in a concrete example in the end. They 
are the lowest dimensional operators that allow RHNs to 
couple to baryon number directly [ 20 ]. 

In order to have successful dexiogenesis, the coefficients 
of the relevant dim -6 operators cannot be too small. To 
avoid excessive low energy baryon number violation for 
A > 1 TeV, those operators would mainly involve right- 
handed third generation quarks, which help avoid severe 
constraints from nucleon decay data, discussed in more 
detail below. To see this, note that quark mass diagonal- 
ization can induce operators that involve light quarks, 
in the presence of left-handed fields. For this reason, 
we will not further consider tree-level dim -6 operators 
NdQQ in Eq. ([2]) . These operators can still be generated 
from Ntbb via radiative corrections, but would not lead 
to severe constraints. 

In light of the above discussion, throughout this work, 
we will focus on operators involving right-handed third 
generation top and bottom quarks (originating from the 
first term in Eq. ([2]), with explicit spinor contractions) 


/^3R 

^BV 


= A^ 


[N c a P R b)[i c P R b] 
A 2 


( 3 ) 


where Pr = (l+ 7 5 )/2 is the right-handed projector. The 
dominance of the third generation could be expected from 
a connection to UV flavor dynamics. Operators with 
other combinations of chiralities and flavors can in prin¬ 
ciple be present, but they must be more suppressed [23] . 

Remarkably, the addition of the operators in Eq. © 
provides all the necessary ingredients encoded in 
Sakharov’s conditions [15] for baryogenesis: (i) these in¬ 
teractions are manifestly baryon number violating, (ii) 
their complex coefficients provide a source of CP vio¬ 
lation, and (in) if the Universe has a low reheat tem¬ 
perature T rh < M a , then the N a , assumed to be non- 
thermally produced throughout this work, will decay out 
of equilibrium. This mechanism, dexiogenesis, allows 
Trh 100 GeV, since the baryon asymmetry is directly 
generated and hence electroweak sphalerons do not need 
to be active. 

Let Ni be the lighter of the two RHNs in our setup. 
Then, the interference of the tree and the 2-loop dia- 


with F(x) = 1 — 8x — 12x 2 logx + 8 a ; 3 — a; 4 . 

In the presence of the higher dimensional operator 
Eq. @ with a TeV scale cutoff, Ni decays induced 
by neutrino Yukawa interactions [Eq. Q] are subdom¬ 
inant, for values of Mi near the weak scale. Given a 
realistic seesaw mechanism for the SM active neutrino 
masses, in general we have y^ < 10 ~ 6 yjM \/(200 GeV) 
in the absence of fine tuning [21] , The induced N\ —> 
Wl decay rate is then estimated to be ~ 

10 ~ 12 GeV (y^-/10 ~ 6 ) 2 [Afi/(200 GeV)]. We find that, 
for Mi of a few hundred GeV and A/v^Ai < 25 TeV, 
that rate is smaller than the baryonic decay rate. 

The baryon asymmetry can be conveniently obtained 
using the unitarity cut method [25] 

Im(A 2 Af) /MA 4 MiM 2 
£ 3072tt 3 |Ai| 2 ^ A ) (M| - M 2 ) ' 1 ’ 

The relation between the above asymmetry and the 
baryon number to entropy ratio rj = n R /s ~ 10 -io m 
depends on the non-thermal production mechanism for 
Ni, but it can plausibly be 77 ~ e/100. For example, let us 
assume that a heavy modulus, such as an inflaton, decays 
equally into radiation and /Vi, which promptly decays. 
We will take the reheat temperature to be Trh ~ 1 GeV. 
Then, one can estimate r) ~ e/g* where g* ~ 100 is the 
number of relativistic degrees of freedom at Trh- Al¬ 
ternatively, if the modulus decays exclusively into Ni, 
and it is the decay of the Ni population that reheats 
the Universe, we end up with rj ~ eT RR /Mi, which for 
Mi ~ 100 GeV, again yields 77 ~ e/100. Hence, for 
Mi ~ M 2 and A a ~ 1, we typically require Mi /A > 0.1. 
Consequently, for Mi < 1 TeV, relevant for collider phe¬ 
nomenology, the cutoff scale must be sufficiently low, 
A < 10 TeV. Let us then examine the experimental con¬ 
straints on A. 

It turns out that nucleon decay limits provide the most 
stringent lower bound on A in the above model where 
RHNs violate both lepton and baryon numbers. While 



FIG. 2: One of the leading diagrams that yield proton decay. 
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the operators in Eq. ([3| do not contain light quarks, 
quantum loop corrections can induce nucleon decay via 
these baryon number violating interactions. Fig. [^pro¬ 
vides a sample two-loop diagram for proton decay, with 
a rate given by 


E(p —► 7 t + v) 


(1 


g A ) 2 a 2 m p _ 2 


32^/2 


( 6 ) 


where gA = 1-27 is the nucleon axial charge, f n = 
131 MeV is pion decay constant, lattice calculations [26] 
yield the form factor a « —0.01125 GeV 3 , and 


^ A q cdG 2 F m,mlV^ d V: b V: b 

K ~ (16tt 2 ) 2 A 2 


( 7 ) 


is the Wilson coefficient from estimating the two-loop 
diagram in Fig. [2] The angle 6 a is the mixing be¬ 
tween N a and the SM active neutrinos. The hadronic 
mass scale A qc d ss 200 MeV must be introduced un¬ 
der a symmetry argument. The operator we started 
with is [N^P R b][t c P R b\, and after the W-loop dressing 
as in Fig. [2j the operator for proton decay turns out to 
be [N^P R d][u c P R d] (which is the radiatively generated 
NdQQ operator mentioned earlier). The fact that one of 
the down quark is still right-handed implies an external 
(constituent) quark mass insertion ~ A qc d- 
The resulting proton decay life time is 


r(p —> tt + is) 


2.5x 10 32 yr 


A/ y/A a 
1.5 TeV 


4 


9a 

10^ 6 


• ( 8 ) 


The current experimental lower limit on the p —>■ 7r + i' de¬ 
cay channel is 1.6 x 10 31 yr [27]. Hence, the above lifetime 
([8]) is not far from the current limit and, in the region of 
parameters considered in our work, can be within the 
reach of future nucleon decay experiments ( 28 ; : 25 ] , 

Here, we also address a potential bound from requiring 
that the asymmetry in Eq. ([5]) is not washed out after 
baryogenesis. The low reheat temperature assumption 
mentioned before can ensure that processes mediated by 
the operators in Eq. such as bb —> Nt , are effectively 

turned off. However, loop processes similar to those de¬ 
picted in Fig. [2] can lead to baryon number violation me¬ 
diated by lighter states. Let us assume, for illustrative 
purposes, that the reheat temperature is Xrh ~ 1 GeV, 
well above the temperature at the onset of Big Bang 
Nucleosynthesis. We then have to make sure that the 
analog of neutron-anti-neutron oscillation mediated by 
(css) (css)/A' 5 is not active. A straightforward compari¬ 
son with Eq. ([7]) yields 


A' 5 


( Aged VtdYub 

U rn s v? s v; b 


2 

Mi, 


( 9 ) 


where m s ss 100 MeV is the strange quark mass. For 
Mi = 200 GeV, we find A' ~ 2 x 10 9 GeV. The rate 


of the css •<-> css process is estimated to be Ta,b =2 ~ 
7rh/A ,1(J ~ 10" 93 GeV, for T RH ~ 1 GeV, which is com¬ 
pletely negligible compared to the Hubble rate at this 
temperature, H ~ T^ H /M planck ~ 10 -19 GeV. We also 
note that bounds from neutron-anti-neutron oscillation 
would not constrain our model, since that process in¬ 
volves up and down quarks, for which the corresponding 
suppression scale is larger than the above A' scale, and 
much higher than the current limit [50] . 

An immediate consequence of Eq. © is the possible 
production of same-sign top quarks at the LHC and fu¬ 
ture hadron colliders, due to the Majorana nature of 
RHNs, as shown in the left panel of Fig. [3] In this pro¬ 
cess, the RHN N a and a top quark are first produced, and 
then N a decays into another top quark and two bottom 
quarks. Because it is a Majorana particle, an on-shell N a 
is equally likely to decay into tbb or ibb final states. The 
violation of baryon number is manifested in terms of the 
violation of top quark number (by two units). The sign 
of the top quark can be inferred from its leptonic decays. 
For a RHN with a few hundred GeV mass and the effec¬ 
tive cutoff scale A/y/X „ of a few TeV, we find that the 
cross section for this process can be as large as ~ 0.3 fb 
in the LHC Run-II at 13 TeV. The main background for 
this signal is from ttbb final states with the lepton charge 
from a top quark decay misidentified, which is suppressed 
by the small misidentification rate m ■ In Table [I| we 
list the leading order cross sections of our signal for sev¬ 
eral sample mass values of RHNs. These points have not 
been excluded by the existing LHC data. For example, 
with M a = 200 GeV and A/y/X^ = 1.5 TeV, the cross 
section at 8 TeV is 0.07 fb, which implies only 1-2 events 
given the existing integrated luminosity ~ 27 fb _ , and it 
is further suppressed by the top quark leptonic branching 
ratios. 


u(pp —> tN —> ttbb ) 

M a 

200 GeV 

500 GeV 

800 GeV 

1 TeV 

yfs = 13 TeV 

0.34 fb 

0.16 fb 

8 x 10~ 2 fb 

5 x 10~ 2 fb 


TABLE I: Same-sign top quark production cross section, 
at the 13 TeV LHC, via a Majorana RHN and the con¬ 
tact operators in Eq. ©■ The cutoff scale is fixed to be 
A/y/K = 1.5 TeV. 


Following the same logic as introducing RHNs to make 
the SM renormalizable, we now discuss a UV comple¬ 
tion that generates the effective operator Eq. (J3|. Given 
a TeV scale cutoff, it is possible to directly probe the 
heavy particles in such a model in LHC Run-II and fu¬ 
ture hadron colliders. The model is an extension of the 
SM that contains a color-triplet scalar, T, with quantum 
numbers (3,1,1/3). The corresponding Lagrangian is 

Ajv = f a TN c a P R b + r T* t c P R b + M 2 |T| 2 . (10) 

In fact, this is the simplest model that yields the flavor 
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FIG. 3: Feynman diagrams for same-sign top quark events that can happen at hadron colliders, using the Majorana nature 
of RHNs and the baryonic interactions of Eq. j3j. Left: pp —» tN production via the contact operator Eq. followed by 
the decay N —> tbb. Right: process in the UV complete model, pair production of color-triplet scalars T,T*, and followed by 
T —> Nb, T* —> Nb, and N —> tbb (via virtual T in parentheses). Baryon number and top quark number are broken when both 
RHNs decay into top quarks using a Majorana mass insertion. 


and color structures of the effective operators in Eq. ^ , 
after integrating out the color-triplet scalar T, corre¬ 
sponding to a cutoff A a /A 2 = /o/'/My. The TeV scale 
cutoff as discussed above can be naturally obtained for 
M t ~ TeV and f a f ~ 0(1). We note that in the above 
UV model, it is possible to have baryogenesis through the 
decays of the T particle m- We will not further explore 
such a possibility in this work. 

The introduction of the scalar T could offer richer phe¬ 
nomenology at colliders. If light enough, T, T* can be 
pair produced at hadron colliders. Each triplet will first 
decay into N + b, which is then followed by subsequent 
decay N —> tbb via a virtual T. The above chain of pro¬ 
cesses are represented by the diagram in the right panel 
of Fig. [3] These together result in same-sign top quark 
final states with many 6-jets. In Table |TT| we give the 
leading-order QCD cross section for the T,T* pair pro¬ 
duction at the 13 TeV LHC and a 100 TeV proton-proton 
collider, calculated with MadGraph 1321 . 


a(pp ->• TT*) 

Mt 

1.5 TeV 

2 TeV 

5 TeV 

10 TeV 

y/s = 13 TeV 

0.16 fb 

0.01 fb 



y/s = 100 TeV 

384 fb 

92 fb 

0.54 fb 

4 x 10" 3 fb 


TABLE II: Pair production cross sections of T, T* via strong 
interaction at the 13 and 100 TeV proton-proton colliders. 


Moreover, an additional distinct signal could be dis¬ 
placed vertices from the decay of RHNs, if we take a 
somewhat larger cutoff scale A/y/X^. In fact, we find for 
Mi = 200 GeV and A/yf\[ > 7 TeV, Eq. Q implies 
a displaced decay length ctn ± 100 /im, which would 
be detectable at the LHC [ 35 ] • This could result from 
Eq. (101 for Mp ~ 1 — 2 TeV and f a ~ /' ~ 0.2. Events 
with same-sign tops and displaced vertices would be quite 
striking and hard to miss in collider experiments. Mean¬ 
while, if the corresponding neutrino Yukawa coupling of 
Ni is y A N > 10~ 7 , sufficient to explain the solar neutrino 
mass difference the partial decay rate of Vi -> Iff 
can be as large as order one. The leptonic decays can be 


used to identify N\ as a RH neutrino (see, e.g., [51]). 

We note that the T particle introduced above has the 
same quantum numbers as the color-triplet partner of 
a Higgs doublet in the fundamental 5 h representation 
of the SU(5) GUT [ 55 ]- This gives the motivation to 
consider a more unified framework for our scenario 1551 . 
Our light color-triplet T cannot arise from the same 5 h as 
the SM Higgs, whose Yukawa couplings are inconsistent 


with those in Eq. (10), within our framework. One could 


introduce a new 5# = (T, D) scalar, where D is a Higgs 
doublet with quantum numbers (1,2,—1/2) under the 
SM. The T couplings in Eq. ( p~0| ) can then come from the 
SU(5) gauge invariant Yukawa interactions, 

£su(5) = fa^Ala^H + Mp5' H 5' H 

+ /o 103535 ^ (24 ff + ug UT )/Ag UT , (11) 

where 1 0 are RHNs and singlets under SU(5), the lower 
index 3 means only the third generation fermions are 
involved, and the bar over a representation means the 
complex conjugate of it. The last term contains a 
higher dimensional operator which after GUT symme¬ 
try breaking, (24 H ) = v G vt diag(2/3, 2/3, 2/3,-1,-1), 
projects out the t c PpbT* operator in Eq. (10) with /' = 


5/o^G(/t/(3Agut)- At the same time it forbids the dan¬ 
gerous operator QLT , thus evading the usual doublet- 
triplet splitting problem [ 37 ]. As a consequence of this 
setup, the quark Yukawa interaction QbpD of the SU(2) 
Higgs doublet D is also forbidden; it still possesses the 
neutrino Yukawa interaction from the f a term. Without 
further splitting of the T and D components, the mass 
of D also lies at the TeV scale with a leptophilic nature. 
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